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We investigate the effect of equation of state and 
convection formalism on the sound speed in solar in- 
terior by using the asymptotic inversion technique of 
Christensen-Dalsgaard, Gough & Thompson (1989). 
In this technique the sound speed difference between 
two solar models is estimated from the known fre- 
quency differences using, 
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u>o is the frequency of p-mode in the reference model, 
oj is the observed frequency for the same mode and 
w = u>/(£ + 1/2). The function H\(w) which is 
related to the sound speed difference between the 
model and the Sun, provides another measure of the 
sound speed difference. Similarly, H^uS) which re- 
flects the differences in the surface layers is better 
suited to test different formulations of solar convec- 
tion theory. The sound speed difference between the 
Sun and model can be found using the relation 
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where a — c/r and a s — a(i?©). As an independent 
diagnostic for differences in the surface layers we also 
consider the scaled frequency difference between the 
models and the Sun defined by A nj ; = Q n ^5LO, where 
Qn.i = E ni i/E (w nt i), E n ,i is the mode energy and 
Eo(u)n : i) is the value of E n j for I — interpolated to 
the frequency w n ,z. 



Table 1: Properties of solar models 
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Conv 
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r d /Re 




Ml 


CM 


OPAL 


0.7275 


0.71643 


153.0 


M2 


MLT 


OPAL 


0.7275 


0.71659 


153.1 


M3 


CM 


MHD 


0.7267 


0.71761 


153.2 


M4 


MLT 


MHD 


0.7267 


0.71776 


153.2 



models use OPAL opacities (Rogers & Iglesias 1992) 
and incorporate diffusion of helium and heavy el- 
ements using the hydrogen abundance profile from 
Bahcall and Pinsonneault (1992) and heavy element 
abundance profile from Proffitt (1994). The con- 
vective flux has been calculated using either the 
usual Mixing Length Theory (MLT) or Canuto & 
Mazzitelli (1991) (CM) formulation. The proper- 
ties of these models with heavy element abundance 
.2=0.019, are summarised in Table 1. These mod- 
els use either MHD (Mihalas et al. 1988) or OPAL 
(Rogers et al. 1996) equations of state. 
The error estimate due to the uncertainties in the ob- 
served frequencies is obtained by simulating 30 sets 
of frequencies where random errors with standard 
deviation quoted by the observers are added to the 
model frequencies. The variance in the computed 
values of relative sound speed difference gives the 
error estimate. The frequencies used are from the 
GONG network (Hill et al. 1996) spanning month 7, 
month 8 and those derived from the averaged spectra 
of months 4 to 7. 
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Figure 1: Relative sound speed difference between the 
Sun and models constructed using OPAL (solid line) and 
MHD (dashed line) equations of state. The hump near 
r/i?© = 0.95 is visible for the MHD equation of state 



For this investigation we construct different solar 
models having identical physics except for the equa- 
tion of state and treatment of convection. All the 



From the tests performed on inversion results it is 
clear that the differential asymptotic inversion tech- 
nique is capable of providing reasonable results in the 
range O.2i? < r < 0.95i? Q . Figure 1 shows the rel- 



ative difference in sound speed between the Sun and 
models Ml and M3. It is evident from the figure that 
the sound speed in model Ml using OPAL equation 
of state is closer to that in the Sun as compared to 
model M3 which uses MHD equation of state. Most 
of this difference in the radiative interior has come 
from the difference in depth of the convection zone 
between the two models. We can also notice a hump 
around r = 0.95Rq inside the convection zone in 
MHD model. Thus it appears that the OPAL equa- 
tion of state better describes the solar material which 
also confirms earlier results (Basu & Antia 1995a; 
Dappen 1996). 

The two formalisms for stellar convection were tested 
by performing inversions using models Ml and M2. 
It appears that the model Ml, which uses the CM 
formulation for calculating the convective flux is 
closer to the Sun as compared to the MLT model. 
However, a definite conclusion can not be made on 
the basis of inversion results as most of the difference 
occurs in the surface layers where the inversions are 
not particularly reliable. 
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Figure 2: The function H2(u>) between the Sun and 
models constructed with different treatment of convec- 
tion. The sold fine represents CM model and the dashed 
line MLT model. 

Figure 2 shows H 2 {uj) as a function of uj/2tt for mod- 
els Ml and M2. We find that model Ml using CM 
formalism has a significantly smaller variation than 
model M2 using MLT formulation. In particular, in 
CM model H 2 (uj) is practically flat for u < 3 mHz 
probably indicating that CM models are closer to 
the Sun than models constructed using MLT. The 
scaled frequency difference between the models Ml 
and M2 and the Sun is plotted in Figure 3. It is 
evident that the scaled frequency difference between 
the CM model and the Sun is smaller than that be- 
tween the MLT model and the Sun. Therefore, we 
conclude that the models using CM formulation for 
convection are closer to the Sun and further confirm 
earlier results (Paterno et al. 1993, Basu & Antia 
1994, 1995b). 
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Figure 3: The scaled frequency difference between the 
Sun and (a) model Ml, and (b) model M2. 
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